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Final  Report: 

Synthetic  Aperture  Radar  Signal  Processing  and  Imaging 
Using  High  Performance  Computing 

Mehrdad  Soumekh 
Department  of  Electrical  Engineering 
State  University  of  New  York  at  Buffalo 
Amherst,  New  York  14260 

This  document  describes  the  progress  on  the  work  performed  for  “Synthetic  Aperture 
Radar  Signal  Processing  and  Imaging  Using  High  Performance  Computing,”  under  the 
DURIP  Contract  F49620-99- 1-0140  for  the  Air  Force  Office  of  Scientific  Research,  for  the 

period  ending  on  9/30/2000. 

The  grant  was  provided  by  AFOSR  to  establish  a  SAR  High  Performance  Laboratory. 
With  this  support,  the  following  items  were  purchased  for  the  SAR  HPC  laboratory: 

1.  16  Origin-2000  300  MHz  Processors 

2.  8  GBytes  of  RAM 

3.  108  GBytes  of  Fiber  Channel  Hard  Drive 

4.  OCTANE  Workstation  with  a  300  MHz  RS-12000  Processor,  and  Enhanced  SSE 
Graphics 

During  the  past  two  summers,  the  PI  also  collaborated  with  the  members  m  the  AL 
Force’s  Rome  Laboratory  (Contact:  Mark  Linderman)  and  Wright  Patterson  Laboratory 
(Steve  Worrell,  Ed  Zelnio  and  Mike  Bryant)  on  various  implementation  aspects  of  SAR  sig¬ 
nal  and  image  processing  algorithms  on  Shared  Memory  Processors  (SMP)  and  Distributed 

Memory  (DM)  HPCs. 

During  the  first  part  of  the  summer  of  1999,  the  PI  worked  with  the  group  at  the 
Rome  Lab.  to  work  on  the  parallel  implementation  of  his  SAR  imaging  algorithms  on 
HPCs  A  subaperture  digital  spotlight  SAR  imaging  method  was  developed  that  was 
suitable  for  implementation  in  both  SMP  and  DM  HPCs.  The  algorithm  was  tested  using 
a  wide-bandwidth  and  wide-beamwidth  FOPEN  SAR  (P-3)  database. 

During  the  second  part  of  the  summer  of  1999  and  summer  of  2000,  the  PI  worked 
with  the  group  at  the  Wright  Pattern  Lab.  to  implement  his  high-resolution  SAR  imaging 
(wavefront  reconstruction)  algorithms  on  realistic  X  band  spotlight  SAR  databases.  The 
subaperture  digital  spotlight  SAR  imaging  method  was  successfully  implemented  for  an 
operational  wide  angle  spotlight  SAR  system. 

The  results  of  these  collaborations  are  documented  in  the  attached  three  papers. 
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Digitally-Spotlighted  Subaperture  SAR  Image  Formation  Using 
High  Performance  Computing 

Mehrdad  Soumekh  and  Gemot  Gunther 
Department  of  Electrical  Engineering,  201  BeU  Hall 
State  University  of  New  York  at  Buffalo 
Amherst,  New  York  14260 

Email:  msoum@eng.buffalo.edu 

Mark  Linderman  and  Ralph  Kohler 
United  States  Air  Force  Research  Laboratory 
Radar  Signal  Processing  Branch 
Rome,  NY  13441 

abstract 

This  paper  is  concerned 

synthetic  aperture.  This  is  achieved  in  conjunct, on  w  h  on  SMP-HPCs  and 

idieital  spotlighting).  Issues  that  are  associated  with  implementing  in  & 

DMP  HPCs  are  discussed.  The  results  using  the  FOPEN  P-3  SAR  data  are  provtded. 

Keywords:  Synthetc  Aperture  Radar.  High  Performance  Computing 

1.  introduction 

This  paper  presents  a  practical  method  for  implementing  the  high-resolution 
SARi^  algorithm.  Jwn  -  the  SAR 

S  .  r  rrsitativplv  lartre  databases  For  real-time  processing  in  an  operational  SAR  system, 
T^temem  puts  restrictions  on  the  minimum  size  of  the  Random  Access  Memory  (RAM)  of  the  computer 
used  for  the  SAR  wavefront  image  formation. 

Meanwhile,  due  to  the  limited  area  and  cooling  restrictions  on  a  radar-carrying  mrcraft,  the  DM-H  s 
_p  the  practical  choice  for  on  board  processing  of  SAR  data.  Unfortunately,  in  a  DM-HPC,  unlike  a  shar 

PHPP  the  RAM  associated  with  each  processor  cannot  be  addressed  by  the  other  processors.  In 
ZZ  the  RAM  available  to  each  DM  processor  (currently  around  64  to  128  MBytes)  »  ■«#« 
o  egpocessmg  some  of  the  recently-acquired  high-resolu.ion  X  band  or  VHF  band  (FOPEN)  SAR 

- The  work  Of  Mehrdad  So'umekh  was  supported  by  the  Office  of  Naval  Research  under  grants r  N00014- 

96-1-0586  and  N00014-97-1-0966,  and  the  Air  Force  Office  Of  Scientific  Research  under  grants  F49620-99- 
1-0140  and  F49620-93-C-0063. 


L  .  ,  .  „„  t.  available  RAM.  Provided  that  the  subaperture  imaging  algorithm  preserves  the 

coherent^brformation  among  the  lower-resolutiou  images,  the  user  could  coherently  add  the  lower-resolution 

images  to  form  the  desired  high-resolution  SAR  image. 

The  challenge  for  this  processing  is  i)  to  extract  the  SAR  signature  of  a  grven  subpatch  from  the 
,  car  Hat*  a  nrocess  that  we  refer  to  as  digital-spotlighting;  and  li)  appropriate  digital  signal 
— gSof  suCrsTa  m:,  yields  "calibrated”  lower-resolution  SAR  ^  ^  - £ 

combined  coherently.  In  this  presentation,  we  address  these  issued  via 
a-  f  fVro  cab  cicrnal  hat  is  ^ased  on  Gabor’s  wavefront  reconstruction  theory,  in  section  A  \\e  outline 

reconstruction  for  stripmap/spotligh.  SAR  systems.  In  Section 
we  show  how  the  lower-resolution  SAR  images  that  are  formed  with  subaperture  data  could  be  coherently 
,  .  ,  ,  th  f  ii  reSolution  image  This  discussion  also  shows  that  subaperture  processing  is  more 

le  arrays  that  are  processed  than  the  full-aperture  processing,  results 
using  a  realistic  stripmap  SAR  database  (P-3  data)  will  be  provided.  Finally,  we  examine  issues  that  ar 
associated  with  the  implementation  of  the  imaging  algorithm  on  a  high  performance  computer  m  Sec  , 

4,  and  present  the  results  of  a  study  that  we  performed  with  an  HPC. 

2.  Subaperture  Digital  Spotlighting 

We  consider  the  problem  of  imaging  a  specific  target  area  in  the  two-dimensionalspatialdomairi 
range  x  and  cross-range  y  from  a  subset  of  a  synthetic  aperture  (subaperture)  within  which  the  target  ar 
is  observable  to  the  radar;  the  synthetic  aperture  could  be  a  portion  of  a  full  aperture  of  a  spotlight  SA 
system,  or  a  subset  of  the  beamwidth  (subaperture)  in  a  stripmap  SAR  system  that  contains  the  desire 

target  area. 

Suppose  the  length,  of  the  subaperture  is  U;  that  is.  «  €  {-L,  I],  where  «  is  the  syn.he.rc  aperture 
domain.  Let  the  coordinates  of  the  center  point  of  the  target  area  with  respect  to  the  center  porn t  of  the 
subaperture  (that  is,  rt  =  0)  be  <X„YC).  The  user  is  interested  in  .magrng  an  area  in  the  spatial  dom 

that  is  spe  An,  A'.  +  An]  and  y  €  [Yc  -  Yo,Y,  +  Yo], 

with  respect  to  the  c-., ter  of  the  subaperture,  that  is,  „  =  0.  Our  task  is  to  develop  a  multi-dimensional 
digital  signal  processing  method  to  extmcl  the  SAR  signature  of  the  above-mentioned  target  area  from 
spotlight /stripmap  measurements,  a  process  that  we  refer  to  as  digital  spotlighting. 

Let  slu.  u)  be  the  measured  SAR  signal  in  the  fast-time  frequency  w  domain  and  the  synthetic  aperture 
u  domain.  To  formulate  the  digital  spotlighting  concept,  consider  the  target  function  which  wforme  via 
narrow-bandwidth  and  narrow-beamwidth  polar  format  processing  at  near  broadside  [S99,  Section  .  j, 

that  is. 

F{kx,ky)  =  sc(*j,u), 


where 


kx  =  2k  cos  0C  and  ky  =  2k,  sin  6C  -  2 kt 


cos2  6, 


In  the  above,  Rc  =  \f~Xl  +  Tc2  and  6C  -  arctan  and 

sc(u>, u)  =  s(*j, u )  exp  [\/ A2  —  JYC  —  u)2] i 

is  the  slow-time  compressed  SAR  signal.  In  this  case,  the  target  function  f(x,y)  is  simply  the  two- 
dimensional  Fourier  transform  of  the  slow-time  compressed  signal,  that  is,  the  slow-time  compressed  SAR 

signal  in  the  (t.ku)  domain. 


U  [S99,  Section  4.1.1,  K  -  - 

^""(xTvn)  or  .he  the  point  spatial  coordinates  [M^J.  ■»•“  “ 

tna- —  and  fcun  «  2fcc  sin  [MO)  -  *c] 

,  c  »d  :macre  that  is,  the  two-dimensional  Fourier  transform  ^(u-.v)  [  sd^-. ”)  ]  ‘ 
in  the  polar  format  processed  SAR  image,  (The  Fourier  transform  with  respect  to 

which  we  refer  to  as  the  polar  format  processed  recor  .tract, on.  f 

the  variable  w  is  an  inverse  one.)  translated  into  the  following 

Note  that  the  above  transformation  from  (r „A)  to 
rectilinear  coordinates  in  the  spatial  domain. 

Xn  =  —  cos (*n  +  Be)  and  yn  =  Mn(*n  + 

here  *.  —  (k)  -  ^  Ug  St  ^ ^ 

fhtgh  a  ^d— al  at  in  the  («,  *.)  domain  with  the  Mowing  passba„d: 

I C-1  cos(d.  +  «t)-Xcl  <  *>  and  |y  si»W+ *c)  -  V,|  <  Vo 

.  -  arcsin  ,  2  is  ,he  nnpnlnr  slow-time  Doppler  domain  a.  the  carrier  frequency  [S99-  Chapter  2]^ 

where  <p  =  arcsin  V2fcc '  lb  y  nnr»nW  d>  domain  is  suitable  for  the 

We  should  point  out  that  this  definition  of  the  ang  ar  s  ow-  tlight  filtering  by  converting 

narrow-bandwidth  SAR  systems.  One  couW  '”Prov  1)  is  the  angular  slow-time  Doppler 

the  SAR  signal  Sc{ u>,  *«)  into  the  (w,*)  omam  w  ^  ^  at  a  fixed  u  or  fc,  one  must  interpolate 

domain  which  varies  with  the  fast-time  freque  y.  d^tabase  is  converted  into  the  (t.o)  domain 

Sc(u,ku)  from  the  ku  domain  to  the*  omain^  ^  ^  ital  spotlight  window  is  then  applied  to  this 
via  an  inverse  Fourier  transform  wi  re®p  data  are  then  brought  back  to  the  (w.*)  domain  via  a 

(<»  domain  database.  The  digital-spothg  ^  ^  }  domain  data  are  interpolated  to  recover 

forward  Fourier  transform  with  respect  tot.  The  re  ffJrcsin  (faA  or  fcu  =  2fc  sin^. 

the  (w,h.)dom»m  data  „„e  needs  to  implemept  one-dioepsio.aHn.erpolat.on 

Thus,  for  the  w.de- bandwidth  SAR J  ^  from  ,he  ^  domain  t0  the  k .  domain.  It  turns  out 

from  the  k,  domain  to  the  f  domain,  an  .  ( for  digital  spotlighting  results  in 

that  using  the  narrow-bandwidth  assumption,  that  s'/  “  are  negligible  effects  especially 

slight  degradations  near  the  edgesoftlitreconsnic  e  ^  a„d/or  completely  remove  these 

^ ^ 

use  the  narrow-bandwidth  assumption  for  dlgita  spot  ig  if  The  reference  S  A.R  signal  is 

The  procedure  for  digital  spotlighting  can  be  summarised  as  Mows.  The  reference  SAR 
de#nC<i  V‘a  s0(n,w)  =  exp[-j2h  +  *Jl 

the  addition  of  the  phase  term  2h  Rc  to  the  reference  sign^  ensures  tM  the  reference  fas, -..me  pom,  Tc 
L  unchanged.  The  digital  spotlight  filter  in  the  («,*.)  domain  ,s  dehned  via 

f  l  for  |%1  cos(<#>  +  «,)  -  AM  <  Xo  and  '  Ytl  <  ^ 

Wd(t,ku)  —  |q,  otherwise. 


Polar  format  processed  reconstruction  with  the  digital  spotlight  filter  is  ' 

/<(«,  fc»)  =  Wd(t,  ku)  [  sc(u,  u )  ] . 

Then,  the  digital-spotlighted  slow-time  compressed  SAR  signal  in  the  («.*.)  domain  is  obtained  by 

Scd(u,  ku)  =  T(t)  [  fd{t,  ku)  | . 

*  M  this  point,  if  the  - 

is,  sample  spacing  conversion  fr  uc  «*)>  Chapter  21.  After  the  optional  zero-padding 

^ J-  * .» * — - 


Scd(«,«)  =  4.)  ]' 


Fin*  the  digital-spotlighted  SAR  signal  in  the  (in, «)  domain  is  constructed  from  the  following  slow-, hue 

dllprtsion  (mixing  with  the  reference  SAR  signal): 

P  ,  .  f  TTUF.Kn.nd  P-3  (FOPEN)  SAR  data.  Figure  la  shows 

In  this  paper,  we  examine  the  processing  ,  j  390  m  (in  range)  and  600 

the  reconstructed  target  spectrin.  for  2,048  slow-time  samplesof  P- <Uta for  a  300  m  (  SJ 
m  (in  cross-, ange),  that  is.  (*.*)  =  (150,300,  m  which  ,  locaud  at  («  - 

suppressLTto  'V'  also  obtain  the 

-  - 

3.  Coherent  Subaperture  Digital  Spotlighting 

The  reconstructed  target  function  in  Figure  lb  is  formed  via  processing  a  subset  of  the  aspect  angle 

their  three  positive  squint  angle  counterparts)  to  obtain  the  full-resolutton  SAR  .mage  of  the  targe  area^ 
For  this  purpose  one  could  combine  the  targe,  area's  spectral  signature  that  appears  ,n  Figure  la  and 
its  five  other  counterparts  (i.e.,  six  spectra,  hands  formed  by  six  different 

an  inverse  two-dimensional  Fourier  transform  to  form  the  full-resolution  spat,^  domain  SAR  ‘mage,  g 
2  shows  the  combined  spectral  distribution  of  the  target  area  from  its  six  subaperture  spec  . 

Note  that  the  slow-.ime  Doppler  bands  (that  is,  the  *,  domain  coverage)  f°r 

st ructions  (e.g„  Figure  la)  are  smalier  than  the  ^ 

fact  the  sizes  of  the  processed  measurement  arrays  and  the  FFTs  that  are  used  to  iorm 
“baperture  images  are  significantly  smalier  than  the  ones  used  for  the  recons —I 

sh”we  consider  the  Nyquist  sampling  rate  for  the  slow-time  domain  which  is  the  smaller,  of  ISM] 

A  •  /I  _  _  1 


A,,  < 


2  (  sin  #max  “  ^min] 


and  < 


2  sin  flmax 

^min 


2  Sin  ^min  ^ 

^max 


where  Amin  and  Ama)(  are,  respectively,  the  smallest  and  largest  wavelengths  in  the  radar  signal,  and 

+  Y0 


are  the  smallest  and  largest  aspect  angles  of  the  target  area  (a  positive  Yc  is  assumed,  and  some  simplifi¬ 
cations  are  used  to  relax  our  notation;  for  the  complete  story,  see  [S99]). 

In  the  case  of  the  P-3  data,  there  are  12,288  (six  times  2,048  samples)  slow-time  samples  spaced 
at  Au  =  .4  m.  Thus,  to  form  the  full-resolution  SAR  image  by  processing  the  entire  12,288  slow-time 
samples,  the  slow-time  Nyquist  criterion  becomes  Au  <  .22  m;  that  is,  the  slow-time  domain  d;  .a  should 
be  upsampled  to  over  22,000  samples  (see  the  discussion  on  PRF  upsampling  via  the  slow-time  compression- 
in  [S99]).  However,  for  the  subaperture  reconstruction,  the  slow- time  Nyquist  criterion  is  Au  <  .92  m  for 
Yc  =  ±409.6  m;  Au  <  .66  m  for  Yc  =  ±1228.8  m;  and  Au  <  .52  m  for  Yc  =  ±2048  m;  In  this  case,  there  is 
no  need  to  upsample  in  the  slow-time  domain.  Iin  fact,  in  our  processing,  we  reduce  the  sampling  rate  in 
the  slow-time  domain  (depending  on  the  value  of  Yc )  since  the  measured  P-3  data  is  with  A„  =  .4  m. 

4.  Implementation  on  a  High  Performance  Computer 

We  have  developed  the  subaperture  digital-spotlighting  and  coherent  imaging  algorithm  using  a  Mat- 
lab  code.  This  code  has  been  successfully  being  used  to  process  the  following  five  different  SAR/ISAR 
databases: 

i.  NRaD's  X-band  ISAR  data 

ii.  SRPs  UHF-band  (FOPEN)  SAR  data 

iii.  NAVAIR/ERIM's  UHF-band  (P-3  FOPEN)  SAR  data 

iv.  ERIM's  X-band  DDB  SAR  data 

v.  ERIM’s  Ultra- SAR  data 

This  Matlab  code  has  been  converted  to  an  efficient  C  code  that  is  suitable  for  implementation  on 
an  HPC.  The  following  is  a  summary  of  the  outcomes  and  observations  in  using  this  C  code  tc  process 
P-3  data  on  a  16-processor  8-GBytes-RAM  SGI  SMP-HPC.  To  improve  the  performance  of  th-  code  f.om 
the  computational  cost  point  of  view,  we  first  considered  its  sequential  optimization;  we  then  studiju  its 
parallel  implementation.  These  are  described  next. 

A.  Sequential  Optimizations 

Fourier  Transforms:  Processing  SAR  data  requires  numerous  one-dimensional  and  two-dimensional 
Fourier  transforms.  Thus,  it  was  important  to  utilize  a  fast  Fourier  transform  code  in  C.  A  great  deal  of 
research  has  been  done  on  this  subject.  We  chose  to  use  the  fftw  library  [Fri].  This  library  determines 
the  optimal  decomposition  of  a  Fourier  transform  on  the  system  in  use  by  performing  a  number  of  sample 
Fourier  transforms,  and  is  easily  parallelizable.  However,  fftw  does  not  automatically  choose  the  fastest 
way  of  implementing  a  transform:  The  user  still  has  to  test  whether  it  is  more  efficient  to  implement  a 
non-contiguous  transform  or  a  transposition  and  a  contiguous  transform,  and  how  many  threads  to  use  for 
a  parallel  FFT.  These  choices  depend  greatly  on  the  size  of  the  data  used  and  will  also  differ  from  system  to 
system.  There  are  some  general  rules  to  follow,  but  no  automated  functions  that  will  aid  in  these  choices. 
Also,  this  library  does  not  implement  an  integrated  FFT-shift  operation,  so  from  case  to  case,  decisions 
had  to  be  made  on  how  best  to  integrate  or  avoid  this  operation.  The  array  sizes  of  the  SAR  data  were 
adjusted  to  allow  particularly  fast  Fourier  transforms  where  possible. 


.  ,YC-L-Y0, 

=  arctan  ( - — - ) 


and  0max  =  arctan  ( 


Yc  +  L 


Combining  Operations:  The  main  focus  in  the  process  of  sequential  optimization  lies  in  the  simplifi¬ 
cation  of  the  computation  and  in  reduction  of  redundancy.  A  clear  indication  that  a  code  is  not  optimized 
yet  is  when  the  program  spends  a  lot  of  its  time  merely  copying  data..  Very  often,  this  can  be  completely 
avoided  by  forcing  in-place  calculations.  Operations  that  inherently  involve  copying,  like  FFT-shift  or 
zero-padding,  can  sometimes  be  combined  with  each  other  and  with  other  calculations  and  thereby  reduce 
the  copying  significantly.  Another  good  example  of  a  reducible  operation  is  performing  several  phase  shifts 
in  one  step,  e.g.,  chirp  signal  calibration  and  motion  compensation.  This  will  significantly  improve  the 
speed  of  the  code  and  is  trivial  from  a  mathematical  point  of  view.  Unfortunately,  it  also  reduces  the 
-  eadabi?  i  ty  of  the  code  because  operations  are  combined  that  are  considered  separate  from  a  radar  imaging 
point  of  view.  In  a  code  that  is  optimized  in  this  way,  it  may  be  rather  difficult  to  track  the  different  parts 
of  each  processing  step  as  they  may  be  integrated  into  different  modules. 

B.  Parallel  implementation 

Two  very  different  ways  of  parallelization  have  to  be  considered:  Parallelism  can  occur  at  the  lowest 
level,  e.g.,  by  using  a  parallel  Fourier  transform,  or  at  the  highest  level  by  running  the  optimized  sequential 
algorithm  in  parallel  for  different  subapertures  or  targets  areas.  Both  of  these  concepts  have  advantages 
and  disadvantages.  Running  several  independent  sequential  algorithms  in  parallel  is  very  straightforward, 
easily  implemented  with  a  script,  and  keeps  all  available  processors  busy.  However,  the  resources  are  not 
optimally  used:  No  use  is  made  of  the  SMP  architecture  because  each  process  works  on  its  own  memory. 
A  small  part  of  the  computation  performed  by  the  sequential  processes  is  redundant  and  could  be  avoided 
if  the  processes  shared  their  results.  No  single  process  can  use  the  large  amount  of  RAM  available  in  a 
SMP-HPC,  each  of  the  p  processes  can  only  use  an  average  of  1/p  of  the  RAM. 

Low  level  parallelism  on  the  other  hand  requires  a  great  amount  of  programming  because  every  single 
part  of  the  program  has  to  be  separately  converted  from  sequential  to  parallel  execution.  The  resulting 
parallel  program  using  p  processors  will  not  be  p  times  faster  than  the  sequential  algorithm:  Creating  and 
synchronizing  threads  involves  a  significant  overhead,  and  not  all  parts  of  the  program  lend  themselves 
equally  well  to  parallelization.  Also,  the  optimal  sequential  algorithm  may  not  be  the  best  algorithm  to 
parallelize.  Other  problems  that  turn  up  are  specific  to  the  target  architecture:  A  single  processor  working 
on  da*j.  in  its  cache  may  finish  its  calculation  faster  than  several  processors  working  on  the  same  data 
and  mutually  invalidating  their  cache  entries.  Because  of  these  difficulties  in  low-level  parallelization,  we 
decided  to  primarily  use  high  level  parallelization  of  a  single  threaded  process  for  imaging  of  large  target 
areas  at  high  resolution.  However,  for  imaging  a  specific  target  area  from  a  specific  point  of  view,  low-level 
parallelization  is  preferable. 

C.  Evaluation 

As  described  in  Section  3,  it  is  possible  to  generate  full  resolution  images  from  lower  resolution  images 
derived  from  subaperture  processing.  We  therefore  have  to  examine  two  parts  of  the  process:  Generating 
the  subaperture  spectral  images  and  combining  these  images  to  a  full  resolution  spatial  image.  The 
parameters  of  our  problem  are  the  target  (sub-)area  size  and  the  size  of  the  subapertures.  The  target 
area’s  dimensions  in  range  and  cross-range  are  2Ao  and  2Yo-  The  given  P3  database  supports  values  for 
A'o  up  to  350m.  It  contains  24  nodes  with  512  slow-time  samples  each.  We  measured  the  CPU  time  and 
memory  requirements  of  the  sequential  algorithm  for  different  target  area  sizes  of  a  broadside  target  and 
different  numbers  of  nodes  in  the  subapertures.  Our  values  for  Xq  were  50,  100,  200,  300,  and  350  meters. 
Our  values  for  Y'o  were  100,  200,  300,  450.  and  600  meters.  We  used  subapertures  containing  1,  2,  3,  4, 
and  6  nodes. 


„  ou  ,  thp  rPIj  time  ner  target  area  used  to  image  that  area.  Figure  3a  shows  the 
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combining  the  subaperture  images  on  the  total  tuning  is  negligi  e. 

The  following  table  shows  the  top  ten  timing  results,  the  parameters  used  and  the  memory  required 
for  each  of  the  two  stages. 


rpT-  number  of  Average  memory  - - - 

i  y  y  nodes  per  subaperture  for  combination 

g)  lu  w  _ 

0.599154  350  600  6  963  1W7 

0.671711  350  600  4  .  656  ^ 

0.722609  300  600  3  839 

0.726007  350  .  4o0  3  ^ 

0.727505  350  600  3  Q69 

0.732468  300  600  6  9  g3g 

0.732560  350  450  4  ®  9g2 

0.748727  300  600  4  613 

0.764735  300  450  3  43  ^ 

0.774671  350  4vJ  6 

an  average  of  500  MB  of  RAM  for  each  process. 

.  ,  y  _  o=niT1  y  _  450m  and  nodes  =  3  yields  the  best  timing  result  complying 

to  the  best  performance  we  could  expect  if  this  part  were  optimally  low-level  parallelized,  allowing 
process  to  use  the  entire  available  memory. 

The  combination  of  the  subaperture  spectral  images  however  is  very  straightforward,  containing  only 
matrix  additions  and  a  2D  Fourier  Transform.  It  is  therefore  easy  to  implement  a  low-  evel  parallelization 
"^k  We  use  all  available  processors  and  are  therefore  no,  limited  to  a  fraction  of  the  available 
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909 


main  memory. 


D.  Performance  on  a  DMP-HCP 

The  implementation  using  high-level  parallelism  for  the  subaperture  imaging  does  not  make  use  of  the 
memory-sharing  capability  of  the  SMP-HPC  other  than  in  passing  the  results  to  the  coherent  subaperture 
combining.  This  kind  of  independent  parallel  execution  of  the  algorithm  could  just  as  easily  be  implemented 
on  a  DMP-HPC.  Notice  that  instead  of  examining  the  average  memory  usage,  we  have  to  consider  the 
maximum  memory  usage  on  a  DMP-HPC. 

The  current  DMP-HPCs  used  on  radar  carrying  aircraft  have  only  64  MB  of  RAM  per  processor.  This 
limits  the  maximum  memory  that  a  process  may  use  for  subaperture  imaging.  From  our  tests,  the  best 
results  under  these  limitations  are  achieved  using  Xo  =  50m,  Yo  =  200m,  and  nodes  =  1,  yielding  a  total 
CPU  time  per  area  of  0.003081sec/m2.  This  it.  about  4.2  times  less  efficient  than  the  implementation  on 
our  SMP-HPC,  therefore  a  DMP-HPC  must  have  4.2  times  as  many  comparable  processors  in  order  to 
process  the  data  in  the  same  time. 

DMP-HPCs  using  128  MB  of  RAM  per  processor  could  process  a  target  area  corresponding  to  Xq  = 
200m  and  Yq  =  200m  with  one  node  per  subaperture,  yielding  a  total  CPU  time  per  area  of  0.001851sec/m2. 
This  is  1.66  times  more  efficient  than  the  implementation  restricted  to  64  MB  of  RAM.  In  the  first  scenario, 
the  coherent  subaperture  combination  requires  only  '48  MB,  it  could  therefore  run  on  a  single  processor. 
In  the  second  case,  188  MB  of  RAM  are  needed.  However,  the  memory  requirements  for  the  coherent 
subaperture  combination  do  not  restrict  the  processing  parameters  because  the  process  and  its  memory 
requirements  can  be  distributed  among  several  processors  with  some  loss  of  efficiency. 

Note  that  the  most  efficient  choice  of  Xo,  Y0,  and  nodes  is  not  necessarily  the  fastest:  If  the  set  of 
input  data  or  the  target  area  of  interest  is  small,  choosing  the  most  efficient  parameters  will  occupy  only 
a  small  number  of  processors.  It  would  then  be  less  efficient  but  faster  to  split  the  problem  into  smaller 
pieces  so  that  no  processors  will  be  idle. 
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ABSTRACT 

v  cap  database  that  was  originally  intended 

This  paper  address  the  problem  of  processing  an  -  ^  we  use  the  approximation-free  SAR 

for  processing  via  a  polar  format  imaging  gon  m.  ^  compensated  phase  history  (polar  format) 

wavefront  reconstruction.  For  this,  the  measure  ssing  algorithm  that  yields  alias-free  slow-time 

data  are  processed  in  a  multi-dimensional  digi  g  P  forrnation  The  X-band  SAR  system  also 

samples.  The  resultant  database  is  used  for  wavehont  SAR  data  are  used  in  a 

X^^—^r^rd  Moving  Target  Indication  (GMTI).  Results  are  provided. 
Keywords:  Synthetc  Aperture  Radar.  Wavefront  Reconstruction,  Monopulse  SAR.  GMTI 

1.  INTRODUCTION 

•WUhtheadve„,o,powerfuUigitaiSi8n.P— 

and  reconstruction  m  Imag'n®  “^  "{e  ,,,  decade,  a  new  SAR  imaging  method.,  known  as  the 

This  is  particularly  true  in  SAR  imaging.  V  -d  high-resolution  and  accurate  coherent 

^etl^ 

^  klhemodernhigh-— 

SAR  wavefront  reconstruct.on.  gAR  g  ems  with  a  relatively  large,  e.g., 

A  prominent  example  is  encountered  in  X  an  P  f  .  sufficient  to  acquire  about  30,000  slow- 
30-degree  look  angle.  Based  on  the  polar  ormat  processin  ,  reconstruction.  the  user  should 

time  samples  for  this  SAR  "d/or  cross-range  samples  for  alias-free  imaging.  ‘ 

Sterns  that  acquire  data  based  on  the  assumptions  of  the  polar  format 

_ _ _ —  , ,  nnnrteA  hv  the  office  of  Naval  Research  under  grants  N00014- 

.S"5 Tporce  Jttce  Of  Science  Research  under  grants 

1-0140  and  F49620-93-C-0063.  im.a<rintr  fCOOl  the  slow-time  samples  are  not 

— for  thc  - 

resolution  image  still  is  about  150,000. 


processing,  provide  a  motion  compensated  phase  history  (MCPH)  data  that  contains  only  30.000  slow-time 
samples. 

In  this  case,  the  challenge  for  the  user  is  to  develop  a  method  for  converting  the  MCPH  SAR  data  to 
alias-free  SAR  data,  and  form  high-resolution  images  from  the  resultant  database  via  the  SAR  wavefront 
reconstruction  algorithm.  This  paper  is  concerned  with  developing  such  an  algorithm,  and  implementing 
it  for  an  operational  spotlight  SAR  system.  The  application  of  the  method  in  constructing  a  statistic  that 
can  be  used  as  a  ground  moving  target  indicator  (GMTI)  in  a  DPCA-based  spotlight  SAR  system  is  shown. 

2.  WAVEFRONT  RECOSNTRUCTION  USING  MCPH  DATA 

A.  Conversion  of  MCPH  Data  into  Phase  History  Data 

Let  s(u,  u)  be  the  measured  spotlight  SAR  signal  (phase  history  data)  in  the  fast-time  frequency  u 
domain  and  the  synthetic  aperture  u  domain.  Suppose  the  spotlighted  target  area  is  identified  via 

x  €  [Xc  -  A'0, Xc  +  Xo]  and  y  £  [Yc  -  Y0,YC  +  Y0], 

where  (Ac,  Yc)  are  the  squint  range  and  cross-range  coordinates  of  the  center  of  the  target  scene,  and  2Ao 
and  2Yo  are  the  size  of  the  spotlighted  area  in  the  range  and  cross-range  domains,  respectively. 

The  MCPH  SAR  data  are  samples  of  the  slow-time  compressed  SAR  signal,  that  is, 

sc(u,  u)  =  «(.*?,  u)  exp  [j2k  %/A"f  -  (Yc  -  u)2]. 

Thus,  the  phase  history. SAR  data,  that  are  samples  of  the  SAR  signal  s(u>,u),  can  be  recovered  from  the 
MCPH  data  via 

s(u,  u)  =  Sc(j;.  u)  exp  [  -  j2k  \fx\  -  {Yc  -  u)2]. 

The  main  point  that  a  user  should  keep  in  mind  is  that  the  slow-time  sample  spacing  for  the  MCPH  data 
(polar  format  processing)  is 

Rc  ^min 
ti  S  777  i 

41  o 

which  is  higher  than  the  slow-time  sample  spacing  of  the  SAR  signal.  Moreover,  due  to  the  imperfections  in 
the  analogly-spotlighted  beam  pattern  (that  is,  the  side  lobes  of  the  radar  radiation  pattern),  there  is  some 
leakage  from  the  targets  that  fall  outside  the  desired  (spotlighted)  target  area  (that  is,  x  £  (Ac— Ao,  Xc+ Ao] 
and  y  £  [Yc  -  lo ,  Yc  +  Y0]);  these  would  introduce  further  slow-time  Doppler  aliasing  in  the  desired  SAR 
phase  history  data. 

B.  Digital  Spotlighting  and  PRF  Upsampling 

Both  of  the  above-mentioned  slow-time  Doppler  aliasing-related  problems  can  be  surmounted  via  the 
digital  spotlighting  and  PRF  upsampling  in  [S99].  For  this  purpose,  suppose  the  length  of  the  subaperture 
is  2 L;  that  is,  u  £  [-L,  L],  where  u  is  the  synthetic  aperture  domain.  Let  the  coordinates  of  the  center  point 
of  the  target  area  with  respect  to  the  center  point  of  the  synthetic  aperture  (that  is,  u  =  0)  be  (AC,YC). 
The  user  is  interested  in  imaging  an  area  in  the  spatial  domain  that  is  specified  via  x  £  [Ac  —  Ao,  Xc  +  Ao] 
and  y  £  [Yc  —  Yq,  Yc  -f  Vo],  with  respect  to  the  center  of  the  aperture,  that  is,  u  =  0.  Our  task  is  to  develop 
a  multi-dimensional  digital  signal  processing  method  to  extract  the  SAR  signature  of  the  above-mentioned 
target  area  from  the  spotlight/stripmap  measurements,  a  process  that  we  refer  to  as  digital  spotlighting. 


This  can  be  achieved  via  the  following  steps.  The  reference  SAR  signal  is  defined  via 


s0(u,lj)  =  exp [—j2k  y/Xj  +  (Yc  -  ti)2  +  j2k  i2c]; 

the  addition  of  the  phase  term  2k  Rc  to  the  reference  signal  ensures  that  the  reference  fast-time  point  T, 
is  unchanged.  The  digital  spotlight  filter  in  the  (t,ku)  domain  is  defined  via  VTd(t,fcu)  =  1  for 

\^cos(4>  +  ec)-Xe\<X0  and  |y  sin(^+  9e)  -  I'd  <  Y0; 

and  zero,  otherwise.  Polar  format  processed  reconstruction  with  the  digital  spotlight  filter  is 

fd(t,ku)  =  Wd(Lku)  JF(Wfll)  [  5c(u;,  u)  j. 

Then,  the  digital-spotlighted  slow-time  compressed  SAR  signal  in  the  (w,  ku)  domain  is  obtained  by 

Scd( V,  ku)  =  F(t)  [  fd(t,  ku)  j  . 


To  upsample  the  SAR  data  in  the  synthetic  aperture  u  domain  (that  is.  sample  spacing  conversion 
from  Auc  to  Au),  an  appropriate  number  of  zeros  should  be  added  to  the  samples  of  Scd(^-ku)  in  the 
slow-time  Doppler  ku  domain  S99,  Chapter  2].  The  Nyquist  sampling  rate  for  the  slow-time  domain  which 
is  the  smallest  of  [S99] 


Au  <  — -  ----- - — r  and  Au  < 

2  ( sin  ^max  -  0min) 


r  2  sin  6n 


2  sin  6„ 
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where  Amin  and  Amax  are,  respectively,  the  smallest  and  largest  wavelengths  in  the  radar  signal,  and 
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=  arctan  ( - - - J 
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are  the  smallest  and  largest  aspect  angles  of  the  target  area  (a  positive  Yc  is  assumed). 

After  the  zero-padding  in  the  ku  domain,  the  digital-spotlighted  slow-time  compressed  signal  in  the 
(uj,  u)  domain  is  formed  via 


scd(u,\  u)  =  T. 


(-1) 

(fc.) 


Finally,  the  digital-spotlighted  SAR  signal  in  the  (u>,  it)  domain  is  constructed  from  the  following  slow-time 
decompression  (mixing  with  the  reference  SAR  signal): 


sd(u,  u)  =  sc<i(a >,  u)  s0(u,u). 


C.  Results 

Figure  la  shows  the  polar  format  processed  image  using  an  X-band  spotlight  SAR  database  (called 
DDB  data);  the  imaged  area  is  about  380  m  in  range  by  380  m  in  cross-range,  and  the  dynamic  range  used 
to  display  the  image  is  from  -60  dB  to  -15  dB.  Figure  lb  is  the  imaged  that  is  formed  via  the  wavefront 


reconstruction  for  approximately  the  same  target  area;  the  dynamic  range  used  to  display  this  image  is  also 
from  -60  dB  to  -15  dB.  The  polar  format  image  in  Figure  la  exhibits  shifting  and  smearing  degradations 
at  the  edges  of  the  380  m  by  380  m  target  area.  In  fact,  the  polar  format  processed  image  begins  to 
show  shift  and  smearing  degradations  for  cross-range  values  that  are  about  80  meters  awav  from  the  scene 
center.  Note  that  for  the  same  dynamic  range  of  -60  dB  to  -15  dB,  the  polar  format  processed  imaee 
shows  stronger  background  clutter.  The  reason  for  this  is  that  the  strongest  reflector  in  the  scene  is  at  the 

range/cross-range  of  (x,  y)  k  (0,  -160)  m  that  is  smeared  and  shifted  due  to  being  relatively  far  from  the 
scene  center. 

It  is  interesting  to  note  that  the  quality  (resolution)  of  the  polar  format  processed  image  at  the  center 
o  the  target  area  (that  is,  (x,p)  m  (Xc.  !'„)]  is  also  worse  than  the  resolution  obtained  via  the  wavefront 
reconstruction  algorithm  in  this  region.  (This  should  not  have  happened;  the  two  methods  should  be 
equally  good  in  the  center  of  the  target  area.)  We  beheve  that  the  polar  format  processor  is  implemented 
with  a  filter  that  suppresses  some  of  the  streaking  and  smearing  effects  of  the  plane  wave  approximation 

at  the  edges  of  the  formed  image.  Unfortunately,  this  filter  also  worsens  the  resolution  of  (contaminates) 
the  polar  format  SAR  image  at  the  center  of  the  target  area. 

Using  the  wavefront  reconstruction  processing  for  the  X-band  DDB  data,  we  are  able  to  form  the 
image  m  Figure  1c  that  is  500  m  in  range  and  1200  m  in  cross-range;  the  dynamic  range  used  to  display 
the  image  in  Figure  1c  is  from  -75  dB  to  -30  dB.  We  should  mentioned  that  the  target  (spot)  size  of  1200 
m  in  cross-range  is  beyond  the  anticipated  3  dB  beamwidth  of  the  radar  antenna.  In  fact,  the  limiting 
factor  for  our  processing  (m  the  cross-range  domain)  is  the  PRF  of  the  radar  system. 

3.  PROCESSING  OF  UNCALIBRATED  MONOPULSE  SAR  IMAGES  FOR  GMTI 
A.  GMTI  Using  Along- Track  Monopulse  SAR 

Along  track  monopulse  SAR  imaging  system  utilizes  two  radars  for  its  data  collection.  One  radar  is 
used  as  a  transmitter  as  well  as  a  monostatic  receiver.  The  other  radar  is  used  only  as  a  bistatic  receiver. 

a?oni’tr7-k  n  i  c  J T  ProcessinZ  alSorith™  of  the  two  monostatic  and  bistatic  databases  of  the 
ng  tra.k  monopulse  SAR  system  to  obtain  two  coherently  identical  SAR  images  of  the  stationary  targets 

in  the  scene.  While  the  stationary  targets  appear  the  same  in  the  monostatic  and  bistatic  SAR  images 
however,  the  same  is  not  true  for  moving  targets.  images, 

tarJ71S/7  iS  7*  ba7  fOF  de7OPing  a  Statistic’  which  we  refer  to  as  the  difference  image,  for  moving 

hi  H  ff  7  ^  m°DOStatic  SAR  imaSe  by  fm(x,y)  and  the  bistatic  image  bv  fb(x,y) 

the  difference  image  for  moving  target  detection  is  defined  via  the  following:  ' 

fd{x,  y)  =  fb(x,  y)  -  fm  ( x ,  y). 

f  (sv) ™TnTv)Z ^  tharthhe  7°  56 °f im7geS  that  ^  f°rmed  in  a  mon°P”b*  SAR  system,  that  is. 
f„(x  y)  and  /6(i  j),  have  a  slight  slow-time  difference.  In  this  case,  the  difference  of  the  two  monopulse 

ges  s  ould  yield  information  for  GMTI  (since  the  stationary  targets  do  not  vary  in  the  slow-time) 
owever.  due  to  miscalibration  between  the  two  radars  of  the  monopulse  SAR  system,  a  better  model  to 
relate  the  monopulse  SAR  images  of  a  stationary  target  area  is 

fb{x'  y)  =  fm(x^  y)  *  *  h(x,  y), 

where  ,,  denotes  two-dimensional  convolution  in  the  spatial  domain,  and  klx,y)  is  an  unknown  (miscal 
i  ration)  impulse  response.  Hence,  one  has  to  perform  a  blind  calibration  of  the  two  images.  A  method 


for  this  using  a  two-dimensional  adaptive  filtering  and  its  implementation  via  a  signal  subspace  processing 
method  are  described  in  [S99,  Chapter  8]. 

B.  Signal  Subspace  Processing 

Adaptive  filtering  methods  have  been  suggested  to  solve  the  above  problem  in  one-dimensional  cases. 
To  apply  these  adaptive  filtering  methods  in  the  two-dimensional  problems,  consider  the  discrete  measured 
data  in  the  domain.  The  impulse  response  h(x,  y)  is  modeled  by  a  finite  two-dimensional  discrete 

filter  hmn;  the  size  of  the  filter,  call  it  (Nx,Ny),  is  chosen  by  the  user  based  on  a  priori  information.  In  the 
following  discussion,  we  choose  both  Nx  and  Ary  to  be  odd  integers,  and  (nx,ny)  —  (Ax/2  —  .5,  Ay/2  —  .5). 
Then,  the  above  continuous  model  is  rewritten  in  the  following  discrete  form 

nx  tty 

=  ^  ^  ^  ^  hmn  fm{xi  Vj  “  TlAy)  +  /e(x,,  J/j), 

m=— nx  n—~ny 


where  ( Ax,  Ay)  represent  the  sensor  sample  spacing  in  the  (x,  y)  domain.  In  the  adaptive  filtering  approach, 
a  solution  for  the  impulse  response  hmn  from  the  knowledge  of  and  fb(xi^yj),  call  it  hmn,  is 

obtained  via  minimizing  the  error  function 

ttt  tty 

I  /*(*«%)-  51  hmn  fm(xi-mAx,yj-nAy)\2. 


The  resultant  solution  is  used  to  estimate  fb(xi,yj )  via 


ttX  ^y 

fb(^i^Vj)  =  y  ^  hmn  fm{xi  “  mAx,  J/j  —  TlAy). 

771—  ~Tljp  71  ••  71  y 

The  statistic  used  for  detecting  the  foregin  object  is  constructed  from  the  following: 

fd(xi,yj )  =  fb{xi,yj )  -  fb(xi,yj). 


In  the  one-dimensional  problems,  the  solution  for  hm  is  formed  via  computing  the  inverse  of  a  large 
covariance  matrix,  a  recursive  LMS  (gradient  descent  adaptive)  algorithm.  These  methods  may  be  utilized 
in  the  two-dimensional  problems  via,  e.g.,  reshaping  the  two-dimensional  arrays  into  one-dimensional  arrays. 
This,  however,  requires  processing  very  large  matrices,  especially  for  the  covariance  matrix  and  the  reshaped 
discrete  filter. 

The  signal  /6(x*,  t/j)  is  the  projection  of  /&(x*,  yj)  into  the  linear  subspace  which  is  defined  by  fm(Xi,yj) 
and  N  -  1,  where  N  =  NxNy,  of  its  shifted  versions;  i.e.,  =  [  fm(xi  —  mAx,yj  —  nAy);  m  = 

— nx, . .  .,nx,  n  =  —ny, . .  .,ny  ].  Thus,  it  is  sufficient  to  identify  the  signal  subspace  $.  and  then  obtain  the 
projection  of  fb{xi,yj)  into  this  signal  subspace  to  construct  fb(xi,yj).  Let  xpe(xi,  yj).  £  =  1,2,  ...,Ar,  be 
a  set  of  orthogonal  basis  functions  which  spans  the  linear  signal  subspace  of  i.e.,  =  [  ipeixi^j);  £  = 

1. 2, . . .,  N  ] ,  where 


<  =  i  o’ 

1  4  '  ^ 


for  £  =  fc; 
otherwise. 


To  generate  this  signal  subspace,  one  can  use  Gram-Schmidt,  modified  Gram-Schmidt.  Householder  or 
Givens  orthogonalization  procedure.  The  size  of  the  signal  subspace,  i.e.,  N ,  depends  on  the  user's  a  priori 
knowledge  of  the  number  of  the  nonzero  coefficients  in  the  discrete  model  of  the  impulse  response  h(x.y). 
For  instance,  if  the  discrete  h(x,y )  contains  (Nx,Ny)  non-zero  pixels,  then  we  should  select  A  =  NrNy. 
In  practice,  the  exact  value  of  NxNy  is  not  known.  In  this  case,  an  estimate  should  be  used  based  on  the 
maximum  anticipated  degree  of  shift  and  calibration  errors  between  the  two  sensors. 

The  projection  of  the  fb{xi,Vj )  into  the  basis  function  J/j),  which  is  identified  by  the  series 

coefficient  at  {l  =  1,2,. . .,  JV),  is  found  via  the  following: 

ae=<  fb  ,  xpi  fb(xi,yj) 

*  j 

The  projection  of  fb(xi,Vj)  into  the  signal  subspace  is 

h{xi,Vj)  =  a(  rj>e(xi,  yj). 

£=l 

The  signal  subspace  difference  image,  i.e.,  the  statistic  for  detecting  the  foreign  object  is  constructed 
via  (6b)  fd(xi,yj )  =  /6(x,-,  yj)  -Jb{ x,*,  yj).  Note  that  both  the  adaptive  filtering  method  and  the  signal 
subspace  projection  of  seek  the  same  minimum  error  energy  solution  for  the  estimate  of  fb{xi,Vj)  in  the 
linear  subspace  of  fm(xn  Vj)  and  its  shifted  versions. 

C.  Results 

We  now  examine  the  GMTI  results  that  are  obtained  using  a  DPCA  (two  antenna)  database  for  the 
X  band  spotlight  SAR  system  of  Figure  1.  Figure  2a  shows  the  SAR  reconstruction  of  a  broadside  target 
area  that  was  obtained  with  the  Channel  1  SAR  data  of  the  DPCA  system.  The  Channel  2  reconstruction 
resembles  the  image  in  Figure  2a  (not  shown  here).  However,  the  coherent  or  noncoherent  (magnitude) 
difference  of  the  images  of  Channels  1  and  2  would  not  result  in  nulling  of  the  stationary  targets  (not 
shown  here).  Figure  2b  shows  the  block-base '^sigrru  subspace  difference  image;  for  this,  we  use  20  pixels 
by  20  pixels  blocks  (approximately  5  m  by  5  n*/-,  and  a  filter  size  of  5  pixels  by  5  pixels.  The  vertical  streaks 
in  Figure  2b  represent  the  signature  of  the  moving  targets.  Note  that  some  of  these  streaks  (i.e.,  moving 
target  signatures)  are  not  visible  in  the  original  SAR  reconstruction  of  Figure  2a. 

Figure  3a  shows  the  SAR  reconstruction  of  an  off-broadside  target  area  that  was  obtained  with  the 
Channel  1  SAR  data  of  the  DPCA  system:  the  target  area  is  composed  of  a  foliage  region  with  moving 
targets  in  its  surroundings.  The  Channel  2  reconstruction  resembles  the  image  in  Figure  3a  (not  shown 
here).  The  coherent  or  noncoherent  (magnitude)  difference  of  the  images  of  Channels  1  and  2  would  not 
result  in  nulling  of  the  foliage  (not  shown  here).  Figure  3b  shows  the  block-based  signal  subspace  difference 
image.  The  vertical  streaks  in  Figure  3b  represent  the  signature  of  the  moving  targets.  Similar  results  for 
another  target  area  that  contains  both  stationary  and  moving  targets  are  shown  in  Figures  4a-b. 
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Figure  lb.  Wavefront  Reconstruction;  dynamic  range  is  -60  dB  to  -15  dB 
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Figure  lc.  Wavefront  Reconstruction;  dynamic  range  is  -75  dB  to  -30  dB 
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ABSTRACT 

This  paper  is  concerned  with  a  representation  of  a  target’s  complex  Synthetic  Aperture  Radar  signature 
that  could  be  used  for  classification  purposes.  In  this  representation,  the  complex  SAR  signature  of  a  desired 
target  area  (chip)  as  a  function  of  the  radar  frequency  and  aspect  angle  (u>,<p)  (i.e.,  the  target  area  phase 
history  data)  are  shown  to  directly  map  into  the  two-dimensional  spectrum  of  the  target’s  image  via  a 
nonlinear  transformation;  the  same  information  base  in  the  (±J,<t>)  domain  is  shown  to  be  retrievable  from 
the  digitally- spotlighted  complex  SAR  signature  of  the  desired  target  region.  For  the  classification  problem 
(Automatic  Target  Recognition),  the  resultant  complex  SAR  signature  in  the  (w.o)  domain  is  compared 
with  a  finite  and  discrete  set  of  reference  complex  SAR  signatures  via  a  process  that  we  refer  to  as  signal 
subspace  matched  filtering. 

Keywords:  Synthetc  Aperture  Radar,  Automatic  Target  Recognition,  Signal  Subspace  Matched  Filtering 

1.  INTRODUCTION 

This  paper  presents  a  new  framework  for  Automatic  Target  Recognition  (ATR)  in  Synthetic  Aperture 
Radar  (SAR)  that  is  based  on  classification  of  the  phase  history  data  of  individual  targets.  In  the  conven¬ 
tional  SAR-ATR  algorithms,  the  magnitude  only  or  the  complex  SAR  image  of  a  target  chip  is  analyzed 
to  determine  the  target  class.  Depending  on  the  relative  coordinates  of  a  target  in  the  imaging  scene  and 
the  flight  path  of  the  radar-carrying  aircraft,  the  SAR  image  (chip)  of  a  target  possesses  a  spatial  warping 
that,  should  be  incorporated  and/or  compensated  for  in  the  SAR  image  formation  and  the  ATR  algorithm. 
Furthermore,  the  formed  SAR  images  exhibit  certain  (slant  plane)  parametric  variations  and  even  aliasing 
errors  and  erroneous  shifts  that  depend  on  the  type  of  the  SAR  imaging  algorithm  used  and  the  individual 
who  developed  it. 

For  the  success  of  a  SAR-ATR  algorithm,  the  spatial  warping  and  slant  plane  parameters  of  the  SAR 
imaging  algorithm  should  be  available  to  the  user.  However,  often  some  or  all  of  the  above  variations 
and  parameters  of  the  formed  SAR  image  are  unknown  to  the  SAR-ATR  user.  These  ambiguities,  in 
conjunction  with  possible  aliasing  and  erroneous  shifts  in  the  formed  image,  result  in  unknown  geometric 
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“STS*  - 

mentioned  problems.  In  the  proposed  approach  S  A1  ^  probability  of  detection  to  identify  the 

algorithm  with  a  high  probab.hty  of  to ^  Next,  the  coordinates  (not  the  ent.re  chip)  of 

coordinates  of  suspected  targets  m  * .  SAR ^imaging  ^  ^  digital_spotiighting  algorithm  that  extracts  the 
each  suspected  target  is  passed  to  *  a  ^  coordinates.  We  wiU  also  show  that  the  SAR 

SAR  phase  history  data  of  a  neigh  or  two-dimensional  Fourier  transform  of  the  desire 

phase  history  data  may  also  be  formed  bv  he  m  ^  two.dimensional  spatial  frequency  domain, 

complex  chip  that  is  Mowed  y  a  non  ine  g.  ^  (phase  history)  of  the  chip  area  at  the  desired 

The  resultant  database  is  the  comp  ex  target’s  complex  radar  signature  with  respect 

coordinates  that  carries  information  on  the  This  database  could  be  passed  through  a  finite 

,0  the  fast-time  (radar)  frequency  “pMtal  *  M  Jgorithm  is  compared  with  a  set  of  preexisting 
order  parametric  modehng  algor, thru.  The  °u,pu t  o  6  the  class  of  the  targe,  a,  the  tuput 

parametric  models  (for  “U  a’suspec.ed  target. 

coordinates  where  the  original  ATR  g  jrfhnr  and  matching  is  introduced,  an  algorithm  that 

A  special  algorithm  for  the  features: 

we  call  siynuf  subspace  matched  fUtenng.  that  possesses  the  following 

i.  The  signal  subspace  matching  algorithm  is  no.  Such' “ S 

(i.e..  variations  of  the  radar  radiation  pa  tern  ,  a  tar-et’s  SAR  signature,  are  one  of  the 

"!g  t  c^fe^R  signature  of  a  target  in  the  Cassihcation  problems. 

.  a.,.-  SAR  signatures  are  commonly  measured  or  generated  (e.g„  ua 

ii.  The  pre-existing  or  reference  arg  e  racit-unie  (radar)  frequencies  and  aspect  angles. 

.  the  X-patch  algorithm)  a,  a  finite  and  d, sere  :  set  c o[  a  urge,  in  the  ,*,*)  domain 

i.e..  («•».  In  this  case,  the  measured  (extracted)  g  gAR  sigriatures.  in  this  case,  the  test 

.  "»•  correspond  to  the  avadabfe  ««« ^ ^  show  Weal  resis(ra,ion.  The  proposed  signal 

and  its  corresponding  reference  SAR  “gna  ”es  misregistration. 

subspace  matching  algorithm  ,s  no,  sens,,,  e  »  U-  type  ^  ^  „  a  function  of  the  radar 

Our  analysis  begins  ”ltk  ,h'  ' -^"rhis  is  followed  bv  identification  of  this  signature  in  the  two- 

-age  (Section  3).  PinaUy,  the  signal  subspace  matched 

filtering  method  is  outlined  in  Section  4. 

"  o™IfsroTTUGHTArN 

PA  PHASE  history  data 

The  type  as  well  as  shape  of  |  “mlmXmLhc  cylinder  is  quite  different  from 

rhosTof  ahtU ThThe  «  site.  ' Thi s  “difference  is  «  bey  feature  mhich  can  be  used  to  d.scnmtna  e  a 

man-made  structure  from  foliage.  .  car  imnnesV  Due  to  the  radar  wave  resonance 

*“ge,s does "ot resemblc 


their  optical  image.  Thus,  shape  (magnitude)  information  in  a  SAR  image  might  not  be  a  reliable  source 
for  target  detection  and  identification.  However,  a  target  cohere*;  SAR  radiation  pattern  or  coherent  . 
image  contains  distinct  features  which  could  be  used  for  target  detect, on  and  identification. 

We  now  present  a  study  to  quantify  the  coherent  SAR  signature  of  targets  A  targe,  exhibits  an 
amplitude  pattern,  which  contains  phase  as  weU  as  magnitude,  when  iUuminated  by  a  radar  system  such 
as  SAR  A  target  amplitude  pattern  ouries  with  the  radar  last-time  .frequency  and  aspect  angle.  For  a 
rada^ located  It  <C,u)  at, a  a  targe,  located  a,  (.„*).  we  denote  the  target  amplitude  pattern  w„h 

an(u.xn,yn  -u).  Note  that  v_u 

arctan  (— - ) 

is  the  target  aspect  angle  with  respect  to  the  radar  broadside. 

The  coherent  SAR  signature  of  the  n-th  target  is  its  contribution  in  the  total  echoed  signal  from  the 
target  scene  When  the  radar  is  located  at  (0.  u)  and  fast-time  frequency  w,  the  coherent  SAR  signature  of 
the  n-th  target,  that  is,  the  tmnsmit-receive  mode  radar-target  radiation  pattern  or  SAR  radiation  pattern 
for  the  n-th  target  can  be  expressed  via  the  following: 

hn{u>.Xn,yn  -  u}=  hr(u^Xn,yn  -  u)  On(w,Zn,Vn  ~  “)  M^n^n  ~  u)t 


Radar-to-Target 


Target 


Target -to-Radar 


in  the  above  model.  MO  >s  the  radar  transmit  mode  radiation  pattern,  MO  the  radar  receive  mode 
radiation  pattern  [S99]. 

.  Substituting  for  MO  and  MO  in  terms  of  their  amPlitude  Patterns  and  spherical  PM  signal  m  t  te 
n-th  target  SAR  radiation  pattern,  we  obtain 

hn(u.  x„,  yn  -  «)  =  aT(w,i„,yn  -  «)  exp [-jky/x^  +  (yn  -  «)'  ] 

.  x  an(ue'.  xn,  yn  ~  ti) 

x  aR(uJ.xn ,  yn  -  u )  exp  [-jky/ x*  +  (j/„  -  u)2  ] 

A  block  diagram  representation  of  the  generation  of  the  SAR  radiation  pattern  for  the  n-th  target  is  shown 
in  Figure  1. 

Using  the  expression  for  the  radar  transmit-receive  mode  radiation  pattern 

h(u>.x,y)  =  Mw-,x,t/)  hR{u;,x,y) 

=  a(~\x,y)  exp(—j2ky/ x2  +  y2) 

where 

a( oj, x,  y)  =  aT(w,  x,  y)  or(w,  x,  y), 
in  the  SAR  radiation  pattern  of  the  n-th  target,  one  obtains 

hn(u;,  Xn,jj/n  -  u)  =  =  °n(w,Xn,!/n  _  u)  xni  Vn  ~  u)  _ _ 

=  a„(ui,  xn,  yn  —  u)  a(u,  xn,  yn  —  u)  exp  [— j2k\Zx^  +  (yn  —  u) 

-  *  N  ✓ 


Amplitude  Pattern 


Spherical  PM  Signal 


We  showed  earlier  that  the  radar  amplitude  pattern  a(u>,xn,yn  -  u)  is  related  to  the  physical  size  and 
to  the  type  of  the  radar  used.  The  magnitude  of  the  radar  amplitude  pattern  I a(u»,z*..y„  -  u)\  dictate^ 


its  power  at  the  coordinates  of  the  n-th  target  at  the  fast-time  frequency  u>.  The  spherical  PM  signal 
eXp  ^  —  j2ky/x^  +  (yn  —  u)2]  is  the  round  trip  phase  delay  of  the  echo  from  the  n-th  target  at  the  fast -time 
frequency  u>. 

Once  a  target  area  (chip)  is  identified  by  a  CFAR  algorithm  as  a  potential  man-made  target,  the  user 
could  extract  the  target  amplitude  pattern  an(-)  via  the  digital  spotlight  algorithm  [S99,  Chapters  5-6];  the 
size  and  location  of  the  digitally-spotlighted  area  are  dictated  by  the  coordinates  and  the  size  of  the  chip. 
This  approach  is  summarized  in  Figure  2. 


3.  REPRESENTATION  IN  SLOW- TIME  ANGULAR  DOPPLER  DOMAIN 
AND  SPECTRAL  DOMAIN-BASED  ATR 

Based  on  the  properties  of  AM-PM  signals  [S99,  Chapter  2],  we  have  the  following  slow-time  Fourier 
transform  for  the  transmit-receive  mode  radar-target  radiation  pattern: 


<u) 


Vn  -  u) 


=  A(u>,ku)  An{ u>,ku)  exp  (  -  jy/Ak2  -  k2u  xn-jku  yn ), 


where  the  target  and  radar  amplitude  patterns  in  the  u  and  ku  domains  are  related  via  the  following: 

an{u,xn,yn  -u)=  An  [2fcsin0„(tt),w] 
a(u>,  xri,  yn  -u)=  A[2fcsin0„(u),u>], 


0n(u)  =  arctan(  — - — ) 

In 

is  the  n-th  target  aspect  angle  when  the  radar  is  located  at  (0,u). 

•  Reminder:  The  amplitude  functions  a(-)  [or  a„(-) ]  and  A(-)  [or  An{-)[  are  scaled  transformations  of 
each  other,  and  not  Fourier  transform  pairs. 

The  angular  slow-time  Doppler  domain  [S99]  :s  defined  via 


k 

<?  =  arcsin(^), 


or 


ku  =  2k  sin  <j>. 

Recall  that  the  slow-time  domain  AM  signal  a„(w,  xn,  yn  -  u )  and  the  slow-time  Doppler  domain  AM  signal 
An(u:.ku)  are  related  via  the  following  scale  transformation: 


An^i^xf) —  „  an(^:.  xn,  yn  u), 

rn 


where  the  mapping  is  defined  via 


ku  =  2k  sin  6n(u). 

Thus,  the  scale  mapping  of  the  AM  signals  in  the  angular  slow-time  Doppler  domain  is 


<j>  =  0n(u). 


r  in  thp  asnect  ancle  On ( tz )  domain,  that  is. 

10 f°™ the am c°mpoMn' of ^  th largct 
,  ,  r.  «  1  Figure  3a  shows  the  aspect  angles  of  the  n-th  target 

This  phenomenon  is  demonstrated m  igure  -  -  >  u  -  v  which  are.  respectively,  6n{-L). 

atn  =  -I.«-  '»* :C1°St  T‘°  Dopplel  mapping  of  the  AM  component  ,n 

£ *«■ ,he slow',ime angular  Dopplc'  * d°ma"'  's 

shown  as  the  polar  angle  domain  of 

kx  =  y/4 fc2  -  k\ 

ky  —  ku 

or  kx  =  2k  cos  4> 

ky  =  2k  sin 

,hich  are  the  range  spatial  frequency  and  cross-range  spatial  frequency  domains,  respectively;  the  above 

is  called  the  SAR  spatial  frequency  mapping  [S99j. 

Next,  we  use  the  slow-time  angular  Doppler  variable  *  to  redefine  the  slow-time  Fourier  transform  o 

the  AM-PM  signal,  that  is, 

S„(w. ku)  =  A„(w.  ku)  exp  (  -  j \M2  -k\x „  -  jku  yn)- 


Via  5„<n>. 2 k  sin o)  =  A„(w. 2 k  sin *)  exp  (  -  j2k  cos 

Also,  we  define  the  AM  signal  of  the  n-th  target  in  the  slow-.ime  angular  Doppler  domain  by 

Apn(~-,d>)  =  A„( w,2fc  sin<p). 

Thus,  we  may  rewrite  the  AM-PM  signal  in  the  angular  slow-time  Coupler  domain  by 

5n(— ',2b  sill  *)  =  At„(w,d)  exp  (  -  ]k,  x^-jk,  »„). 

The  equivalence  of  the  aspect  angle  *„< »)  domain  and  the  slow-time  angular  Doppler r  *  domain  in  the 
mapping  of  the  AM  signature  of  a  targe,  plays  a  key  role  in  a 

crucial  for  target  recognition  or  detection  in  SAR J”  f  ^f^could  tom  its  two-dimensional 
targe,  chip,  call  .t “^dimensTcTal"^^)  Fourier  transform.  Finally,  the  target  complex 

Ukl M  from  the  (b„ by)  domain  into  w.d)  domain  v,a  the 

following  nonlinear  two-dimensional  transformation  (SAR  inverse  mapping). 


o  =  arctan 


(t) 


The  resultant  database  if  then  used  for  classification;  this  procedure  is  outlined  in  Figure  4. 


4.  SIGNAL  SUBSPACE  MATCHED  FILTERING 

,  .  4  ,  thp  reference  and  test  SAR  signatures  of  a  target.  In 

Next,  we  consider  the  problem  of  regis  ra  .  in  the  two  complex  SAR  signatures.  These  subtle 
practice,  the  peer  faces  subtle  phase  and  gate  for  tU  m  SAR  sjrstems.  and  incremental 

variations  are  due  to  the  changes  o  e  P  respect  to  the  radar  (small  shift,  scaling  and  rotation), 

changes  in  the  relative  orientat.on  of  the  target  with  respect 

These  factors  are  outlined  in  the  following. 

•  mnv  varv  slightly  in  the  reference  run  (or 

*■  r  ^  “d  ro,a,‘on  cf  iarsets  in 

,  ri  —  c— „  errors^This  This 

"L  target  SAR 

iV'  finUt^and tdhrcrete  se^oOadar  frequencies  and  aspect  angles  tha,  may  no, 

correspond  to  the  measured  SAR  signature  of  a  test  target. 

,  flKv  4  Ju  6)  and  the  (discrete)  test  target  signa- 
We  denote  the  (discrete!  reference  target  “*"“"**  •„  commonly  available  or  generated 

ture  by  ArM «,*).  (Note  that  a  simulated  migh.  be  measured  only  within 

fo,  all  aspect  angle  J  *  a“0ve-men.ioned  calibration  and  misregistration  errors,  the 

<  Apn 2(W,©)  •  Apnl(w,<p)  >f 

would  not  result  in  ther  ideal  correlation  results.  However, th*  ^  this 

be  circumvented  by  the  X-S UH 

tt«  r  ^  subspace  ^ 

to  be  1  c'ollatii  of  this  subspace  projection  and  the  reference  targe,  s.gnature;  ,.e.. 

<  Apn2 (w,  O)  ,  Apnl(wsd>)  >  • 
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